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Abstract 

In this work, a two-dimensional one-band Hubbard model is investigated within a two-pole 
approximation. The model presents a non-local attractive potential U{U <Q) that allows the study 
of d-wave superconductivity and also includes hopping up to second-nearest-neighbors. The two- 
pole scheme has been proposed to improve the Hubbard-I approximation. The analytical results 
show a more complex form for the gap A(T), when compared to the one obtained in the latter 
approximation. Indeed, new anomalous correlation functions associated with the superconductivity 
are involved in the calculation of A(T). Numerical results in a range of temperatures are presented. 
Moreover, the structure of the quasiparticle bands and the topology of the Fermi surface are studied 
in detail in the normal state. Connections with some experimental results are also included. 
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I. INTRODUCTION 



Superconductivity in strongly correlated systems is a field with plenty of challenging 
problems. Several non-usual properties of high temperature superconductors (HTSC) l| 
still needed to be properly clarified. In particular, some experimental systems exhibit im- 
portant deviations from the standard BCS theory. For instance, the superconductor gap 
behaviour as a function of temperature in some borocarbides display a non-monotonic fea- 
ture at lower temperatures. In fact, such gap decreasing in RNi2B2C{R = Dy, Ho, Er, Tm) 
and ErNi2B2C |2| when the temperature decreases towards to T = 0. This effect ascribed 
due to a competition between the superconductivity and antiferromagnetic correlations, 
which are absent in a weak coupling BCS-like superconductivity, appear naturally from the 
formulation which is shown in this work. Another quite interesting feature of HTSC sys- 
tems is the behaviour found in some cuprates ^ in the low doping region, the so called 
pseudo-gap region, which gives rise to an anomalous Fermi surface leading to a pseudo-gap. 

In this very complex problem, a number of theories has been proposed in order to explain 
the presence of a pseudo-gap region In the present formulation, we claim that the 
appearance of a pseudo-gap can be ascribed to a more detailed many-body treatment in 
which superconductor and AF correlations compete. 

Although BCS-like approach ^ has been widely used to describe these physical systems, 
it is well recognized that superconductivity is a two-dimensional problem in which strong 
correlations play a fundamental role js]]. Thereby, we apply a two-pole approximation js-Sl 
to deal with the strong interaction coupling. Here we consider a d-wave symmetry gap and 
therefore, a non-local attractive interaction is used jo], [lo|. The net attractive interaction 
{U < 0) may result, for example, from the elimination of the electron-phonon like coupling 
through a canonical transformation ll| or, alternatively, from an electronic mechanism 
proposed by Hirsch [l^ which may produce, for a certain range of parameters, an effective 
attractive interaction. 

In this work, we focused on the many-body renormalized normal state of these systems. 

nn 

Our obtained Fermi surface is consistent with recent claims in the literature p|, 113[ about the 
presence of hole-pockets due to antiferromagnetic correlations. Moreover, we discuss some 
thermodynamical properties of the superconducting regime, namely the critical temperature 
Tc, the zero temperature superconducting gap Aq and temperature dependence of the gap 
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A(T) for various dopings S = 1 — ut (with ut = + U-a) and interaction U . The 71^ 
represents the average occupation per site of electrons with spin cr =f , J,. 

This paper is organized as follows. In the section [Tll we present a general formulation 
describing the model as well as the ingredients of the normal state, e.g. the quasi-particle 
and the special characteristics of the Fermi surface (FS). In section llllt we present the 
superconducting state which appears from the application of the two-pole approximation. 
In section [TVt we exhibit self-consistent numerical results and conclusions for both the normal 
and the superconducting states . The Appendix |A] briefly describes the main points involved 
in the two-pole approach whereas in the Appendix [B] the correlation functions involved in 
the Green's functions governing the superconducting and the normal states are displayed. 



II. GENERAL FORMULATION 

The Hamiltonian studied here is, in a standard notation 



where ((...)) indicates the sum over the first and the second- nearest-neighbors of i and /i is 
the chemical potential. The two-dimensional dispersion relation is given by: 

= 2t{cos{kxa) + cos^kyo)) + 4^2 cos{kxa) cos{kya). (2) 

In the present work, we adopted the two-pole approximation which consists in 

choosing a set of operators describing the most important excitations of the system. The 
details of the method are given in Appendix |Al In the present case, the set of operators 
considered is ^^di^a,nf-^di^a,dl_^,nf^^dl_J^. The first two are associa;ec with the normal 
state whereas the last two are associated with the superconductivity 3, Following the 
method exhibited in the Appendix |A1 the one-particle Green's function for the normal state 



is: 



GZik, u) = -^IB^ + ^2fe^ (3) 



with 



lk(T Zka 



1 , U-2U,-Ei + Wi^ 



= 2 + 2xf 

ka 

^2ka = ^Ika (5) 
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and 



The quasiparticle bands are: 



^Ika - 2 2~' 



where 



and 



Xe. = \^(U-ei + WiJ^ + 4(7, (£j - WiJ + U (9) 



g^ ^^f--^-' (10) 

The effective interactions f/i, ?72 and the band shift W^^ are defined in the Appendix \M 
(see Eqs. (IA5P - flA7P ). Here, as we are assuming a paramagnetic state of a translationally 
invariant system, {ni^„) = (ui^^^) = (n^^). It should be noticed that, due to many body 
effects, in the pole structure of the Green's functions, in the normal paramagnetic phase, 
there is a spin-spin correlation function which exhibit only antiferromagnetic (AF) short 
range correlations. That is not in contradiction in our previous paramagnetic assumption. 
Moreover, in order to simplify the notation, we write (n^cr) = n^^r- 

From the Green's function in equation (|3]), we find the spectral function 

A^{k,u) = --lm[GZ{k,u;)]. (11) 

The Fermi surface is obtained from Acr{k,uj = 0). 

III. THE SUPERCONDUCTING STATE 

In the superconducting state, the Green's function G'^'^ is written as: 

Giiik, uj) = ^-L^ (12) 

where 

A'{ijj) = ao + aiuj + 021^^ + a^oj^ (13) 
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with 



ai = 2n^„nEuEi2 - {n + 2n^^)n^^Ef2 - n^-a^u 
-E22[E22 + 2n^^{n^^Eu - 2^12)] 

2 Q ^ ^2 



Q!2 
"3 



?En 



n 



(14) 

(15) 
(16) 
(17) 



and 77, = (1 +n_o-). The Enm are elements of the energy matrix flA2p and the quantity P{uj) 
is defined as: 



P{u) = [{tu - En)in^.uj - E22) - {n^.co - E 



12) 



X [{tu + Eu){n^^uj + E22) - {n^^uj + E12) 
+A|(Ai - A2UJ) 



with 



Ai = aQ + ai— + 



a2 + ^l{l+n\^f 



and 



2 ^ \2 I 2 /I 

r^„^-^) +ri_^(l 



A2 = {1 

The quantities oq, ai, 02 and 6*, are given by 



U 



ao 
02 



^22 



= 2Ei2(2n_,^ii-^i2) 
4n_,Ei2^22 + 2^2^ (^^2 + ^11^22) 



and 



(18) 

(19) 
(20) 

(21) 
(22) 
(23) 

(24) 



e = tnoi.-UiDoi. + 2{S',S'^)) 

with the correlation functions noio-, -Dqio- and {S^S^) defined in the Appendix [Bl 

The main reason that we are adopting the d-wave symmetry is that we are following 

reference 0] where it is claimed that for a large number of HTSC material d-wave gap 

symmetry is the most relevant. 

Moreover, the d-wave symmetry follows also from the fact that in our Hamiltonian we 

consider an attractive delocalized interaction term. Actually the extended s-wave symmetry 

is more favoured for an attractive local interaction as discussed in 9i. 



For the particular case of pairing with rf-wave symmetry, the gap function is 



= 2A[cos{K) - cos{ky)], (25) 

where A is the gap function amphtude. Following the procedure described in references 
71, |8|, the self-consistent gap function has been obtained from the Green's function: 

GfA.) = -M^±ft^£!l (26) 

in which, 

A = -29Aj J2[cos{q^) - cos{qy)]^Fig-^ (27) 

and 



L 



The /So and (ii are 



/9o 



n\{l 



,2U 



/3i = El - [El - 2E^E2 + A|(l 



n 



(29) 
(30) 



with El = Ei2 — n_o--Eii and E2 = E22 — n_„E 
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IV. SELF-CONSISTENT RESULTS AND CONCLUSIONS 

Firstly, we discuss the numerical results for the normal state (T > Tc). Figure [T] shows 
the FS for four different doping values 5. In (a), 5 = 0.30, a well defined electron-like FS 
is found. However, in (b) when 5 is decreased to 5 = 0.20, the topology of the FS changes, 
with the emergence of a hole-pocket enclosing the nodal point (f , f). As a consequence, due 
to low spectral intensity, a pseudogap appears near the antinodal points (tt, 0) and (0,7r). 
Further decrease in 5 intensifies the presence of the pseudo-gap as shown in (c)-(d). This 
unusual behaviour is due to the strong antiferromagnetic correlations coming from the band 
shift Wi^ defined in equation flA7p . 

The effects described above, are corroborated by the quasiparticle band calculation ex- 
hibited in figure [2] which displays the quasiparticle band for distinct doping values 6. In 
fact, while in the higher doping regime the quasiparticle bands cross the Fermi level near 
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(a) nT=0.70 (,,,) (b) nT=0.80 (,,,) 




FIG. 1: (Color)The spectral function A{k,Lj = 0) representing the Fermi surface for different 
dopings 5 = 1 — nx- The model parameters considered here are U = 8t, t = —1.0 eV, t2 = 0.3|t| 
and ksT = 0.1\t\ {ks is the Boltzmann constant). 
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FIG. 2: (Color) The quasiparticle bands intercepted by the chemical potential (// = 0). The model 
parameters and the temperature are the same as in figured) 

(|-, |) and near the antinodal point (0,7r), in the lower doping regime the quasiparticle 
band crosses the Fermi level only near the nodal point (f, f). Such a behavior gives rise 
to a pocket around (f , |). On the other hand, as the quasiparticle band does not touch 
the Fermi level near (0,7r), a pseudogap emerges in that region. As far as we know, the 
emergence of the pseudogap due to an attractive U in a. strong correlation regime, is for the 
first time presented here. The kink observed near the (vr, vr) point of the quasiparticle band 
is caused by the strong antiferromagnetic correlations associated with {Sj ■ Si), which are 
maximum in Q= (vr,7r). The Q is the antiferromagnetic wave- vector. 

Now we discuss some thermodynamical properties associated to the superconducting 



7 



U=8t t=-1.0eV t2=0.3|t| 



1 


' nj=0.90 — 

nj=0.80 




0.669 


1 1 1 1 
'^m nj=0.90 








0.666 










0.663 


1 1 l" 1 








0.528 


1 1 1 1 

nj=0.80 


1 


1 1 \ 


t 


0.525 
0.522 


% 

1 1 1 h 



1 2 3 4 5 0.2 0.4 0.6 0.8 

kaT/ltl keT/ltl 

FIG. 3: (Color)The main figure shows the gap function ampHtude versus the temperature for 
U = 8t, t = —1.0 eV, t2 = 0.3|i| and two different occupations ut- The small figures show the 
regions of low temperatures where the gap presents an unusual behavior. 

state. In figure [3] we describe the gap function amplitude A(T) for U = 8t and two different 
occupations ut, in the lower doping regime. One sees that for a given U (in a character- 
istic strong coupling regime | j| >> 1), the zero temperature gap for ut = 0.90 is higher 
than the corresponding one for Ut = 0.80. Furthermore, the temperature where a non- 
superconducting phase arises is higher for tit = 0.90, i.e., TdriT = 0.90) > TdjiT = 0.80). 
It should be noted that in both cases for riy, in the region of low T, there is a increase in 
the value of the gap amplitude as compared to the zero gap amplitude value. This unusual 
behavior is due to the effect of the strong correlations, since that in the BCS weak correlated 
regime, A(0) is always greater than A(T). 

We stress that in our case this non-monotonic behavior at low temperatures is mainly 
due to correlation functions in the pole structure of the superconduction Green's function 
(see equations ( fT2|) and (l26ll ). To be more precise, we have found in our self-consistent 
calculation a complex interplay between the SC gap behaviour and the AF type short range 
correlations. 

Figure m displays the value of the gap amplitude A(T) as a function of temperature for 
several values of t/, in the strong correlation regime (|y| >> 1). In all cases, for ut = 0.90, 
we note that when U increases, Tc increases also, and the same unusual behavior for A(T) 
appears for very low T. We have calculated the gap function amplitude for several values of 
f/, for different ut {ut = 0.80 and 0.70), and the same behavior is observed. 

A very interesting result is shown in figure [51 Here, we plot ^max - Aq for several values 
of One noted that when U increases, ^max - Aq increases also, until a special value (in 
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FIG. 4: (Color)The main figure shows the gap function amplitude versus the temperature for 
riT = 0.90 and different values of U. The remaining parameters are identical to figure El The small 
figures show the regions of low temperatures for two values of U. For low temperatures, the gap 
presents an unusual behavior. 
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FIG. 5: (Color)The difference between A^ax and Aq as a function of 1^1- The inset shows as 
and Aq have been obtained. 



our calculation I j | ~ 



8) where A^ax - Aq is a maximum and then, A^ax - tends toward 
zero in a very high value of U. When such behavior appears, one has attained the \U\ oo 
limit. Moreover, it should be noted that this unusual behavior occurs for a critical low value 
of \U\, (in our calculation | j| ~ 4), which is a signature of the onset of a characteristic 



1 1 



>> 1, signaling the appearance of a strong correlation regime. 
Our results are in qualitatively agreement concerning the Fermi surface in underdoping 
region with other approaches using the t-J model jl^]. The reason of such agreement is that 
in both approaches the spin spin correlation functions renormalize the band structure giving 
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rise to the appearance of hole-pockets. 

Finally, in order to complement our present calculations, we need to discuss the higher 
doping regime nr ^ 0.80. Moreover, a detailed discussion of the effect of external pressure, 
which affects mainly the ratio y 15|], is needed. These further calculations, are now in 
progress. 



Appendix A: Two-pole approximation 

In the present two-pole approximation js, 7|, the Green's functions are defined as: 

G{uj) = N(a;N- E)-^N 



(Al) 



(A2) 



where, E and N are the energy and the normalization matrices given by 
and 

Nrrm={[An,Al]i+)). (A3) 

In equations f lA2|) - (lA3p . [...,...](+)_ denote the (anti)commutator, and (...), the thermal 
average. The set of operators {An} must satisfy, within some approximation, the relation 

For the normal state of the model ([T]), the energy matrix is given by. 



E 



with 



- /i + f/2 £k^l - /i + f/2 + nW^:^ 

[/i = 2f/5^(nz,_.), 
?72 = 2f/5^(nz,_,n,,_,) 



(A4) 



(A5) 
(A6) 



and n = n_cr(l — ri-a). The correlation function Dn^^ = {rii^^ni^^) is defined in equation 



f lBSp . The band shift Wj^^, can be written as: 



nWf 



ka 



■JO- 



(A7) 



with riijfj, rriija and hj^, defined below. 
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Appendix B: Correlation functions 

The correlation function riija = {dl^dja), is given by: 

'''''' = ^ 5 / ^)duJ (Bl) 

k 

with G^*^ defined in equation (I12p . Assuming i = and toj = t for the 2; nearest-neighbors, 
only one value of riojo-, namely riQio-, is necessary. Considering the same for the second- 
nearest-neighbors t2, we have: 

By using the original Roth's scheme j^, the correlation function Dij^r = {nifjUj^), is 
calculated and written as: 



Dija = n, ' _n (B3) 

Pa, a Pa,— a 

with mij^ = {dl^rij^^dj^) given by 

1_ ^ y e^^^-(«^-^')/(a;)G-/(^, (B4) 



k 

and 



Tlijcr Toiler /T->r\ 

ftija = -i — (B5) 

1 - 

1 - n„„ 

The -Dqio- can be obtained assuming again 2 = and toj = ^ for the 2; nearest-neighbors, as 
it has been done in rioio-- 

The Green's function Gff in IH is: 

G7.\k. ^) = — m 

where 

= 70 + 71^^ + 72^^ + 73^^ (B8) 
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with 



70 = (^?2 - ^11^22) [^3 + n^a{Ei2 - Eu)] 

71 = n^„Eu[Ei2{l + Sn^a) - n^„{Eu + nE22) 

+E22E3 + n„,Ei2(3^3 - n-^Eu) 

72 = n^^{l - n^„fEi2 
73= n^^(l-n_^)2. 



The quantities E^ and n are: 



Ez = E22 - E12 and n = 1 + n_cr. 



The term A'" introduced in equation (1B7P is defined as: 



A"'{uj) = io{l + ni-) + En + 7rN_.(Ei2 + E 



11) 



E- 
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9' '^ 9 

The denominator of the Green's function G"^'^, is given in equation ( IT8|l . 
The term /ijo- presented in the band shift flA7l) . is given by: 



^j(T — -Bjo- + {Sj ■ Sq) 



with 



B 



(SjSq) 



1-/3. 



and 



1 



(5, -80) = - {{SIS,) + (SrSl)) + (5|5o^). 



(B9) 

(BIO) 

(Bll) 
(B12) 

(B13) 

(B14) 



(B15) 



(B16) 



^,^u.. ■ .-^.-^u,- (BIT) 
Particularly, in the paramagnetic state, (S^Sq) = (SJ'Sq), then, (5^ • Sq) can be written 



as: 



{Sj ■ Sq) — {S^Sq ) + (SjSq) 



where, 



and 



{S^Sq ) — {d}j^dj^adl_„dQa) 



1 + 



(1 



d\2 ^jcr'^OjcT ~l~ f^ja^jij) 



(B18) 



(B19) 



(B20) 
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The correlation functions n, m, a and /3, are defined in equations (IBip . (lB4p . (IBSP and (lB6p . 
respectively. 

In order to calculate rijjo- and mjjo- in the normal state, it is necessary to consider A = 
in the Green's functions Gf^ and G^^^ defined in equations f[T^ and flB7l) . In this case, 
Gf. ^ G*^, and G-/ ^ G-^. 
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